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A variety of bursting and spiking patterns arise in mod- 
els for respiratory neurons [1] and other neurons asso- 
ciated with hormone release [2]. These models often 
feature quantities evolving on distinct time scales, such 
as fast voltage and slower ion current activation or inac- 
tivation variables. Furthermore, such systems may com- 
bine multiple interacting oscillatory mechanisms, such 
as intrinsic calcium oscillations together with a calcium- 
dependent, voltage-gated membrane potential oscillation 
mechanism. 

Motivated by the activity observed in these models, the 
goal of this project is to understand bursting dynamics in 
three time scale systems. We are particularly interested 
in knowing how much of the complication in burst pat- 
terns results from the presence of three or more time- 
scales in the system and how much is, rather, a reflection 
of the models being relatively high dimensional nonlinear 
systems with many parameters. With this motivation, we 
construct a model consisting of two copies of Morris- 
Lecar equations with three time scales. By considering 
two viewpoints within the realm of geometric singular 
perturbation theory, we explain the mechanisms underly- 
ing the bursting dynamics of our model system, making 
use of computed critical manifold, superslow manifold 
and bifurcations. To elucidate which characteristics truly 
represent three time scale features, we investigate certain 
reductions to two time scales and the parameter depen- 
dence of solution features in the three time scale 
framework. 

Our results in this area may be useful for characterizing, 
developing models of, and analyzing models of experimen- 
tal data. Moreover, they provide useful information for 
future modeling studies, where a determination needs to 
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be made about how many time scales to include in a 
model to represent some given experimental data. 
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